Spin-orbit-driven BCS-BEC crossover in a narrow Feshbach resonance 
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In this paper we investigate the spin-orbit-driven BCS-BEC crossover physics in a narrow Feshbash 
resonance, and mainly focus on the bound state and the striking molecular signature based on a 
generalized two-channel model. For a weak spin-orbit coupling, only one bound state occurs. More 
importantly, we find two and more bound states in the case of the strong spin-orbit coupling. This 
bound-state number depends strongly on the strength, and especially, the type of spin-orbit coupling. 
The existence of the multi-bound states is a new phenomenon induced by spin-orbit coupling. In 
addition, we find that the molecular number does not vanish (but approaches zero) for a fixed 
positive detuning in the presence of spin-orbit coupling. When varying the strength of spin-orbit 
coupling in the positive detuning, an unconventional molecular signature that the molecular number 
increases and then decreases is revealed. We believe that in experiments our predicted molecular 
behavior is a good signature to detect the spin-orbit-induced superfluid physics. 
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Owing to their unprecedented level of control and pre- 
cision, the investigation of spin-orbit (SO) coupling in the 
neutral atoms has attracted much attentions [lj. Espe- 
cially, a one-dimensional (ID) equal Rashba and Dressel- 
haus SO coupling has been first realized in the ultracold 
87 Rb atoms by a couple of Raman lasers |2|. Recently, 
the same SO coupling has been also created experimen- 
tally in degenerate Fermi gas with 40 K Q and 6 Li In 
the presence of SO coupling, the degenerate Fermi gas 
has exhibited the interesting superfluid physics in both 
three 0-[2l| and lower dimensions [23 - l29j . For exam- 
ple, by increasing the strength of a Rashba SO coupling, 
the Cooper paring gap can be enhanced significantly at- 
tributed to the increased density of state at the Fermi sur- 
face. More importantly, this system may change from the 
Bardeen-Cooper-Schrieffer (BCS) superfluid to the Bose- 
Einstein condensate (BEC) with a new two-body bound 
state called Rashbon [EIH, even on the BCS side with a 
negative scattering. When the effective Zeeman field is 
applied, the 2D degenerate Fermi gas with the Rashba SO 
coupling has an exotic topological superfluid supporting 
the Majorana fermions [23|, which is the heart for realiz- 
ing the topological quantum computing (30| . Recently, a 
universal midgap bound state in topological superfluids, 
which is induced by either nonmagnetic or magnetic im- 
purities in the strong scattering limit, has been predicted 

For illustrating the SO-driven BCS-BEC crossover 
physics, a one-channel model, in which only the atoms 
tuned via Feshbash-resonant technique arc taken into 
account, has been introduced in the previous consider- 
ations |5M29|. However, in degenerate Fermi gas, this 
one-channel model is valid for a broad Feshbash-resonant 
regime T ^> 1 



[32J, where V = ^32 mf i B al g Bl/(h 2 7rE F ) 
is the dimensionless parameter, /ib is the Bohr magne- 
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ton, m is the atom mass, ab g is the background s-wave 
scattering strength, B w is the resonant width and E-p is 
the Fermi energy. In fact, to get a more realistic and 
complete description of the BCS-BEC crossover physics, 
especially in a narrow Feshbash-resonant limit r « 1, 
a two-channel model, which includes both the atoms in 
the open channel and the molecules in the closed chan- 
nel, must be introduced [33M36j |. In the narrow Feshbash- 
resonant regime, some fundamental properties can be ob- 
served experimentally by detecting the striking molecu- 
lar signature [37j , additional to measuring the superfluid 
pairing gap applied usually in the one-channel model [38j . 
More importantl y, n ew quantum phase transitions can 
be predicted (39l.]40j. due to the existence of extra U(l) 
symmetry for the molecular field. 

On the experimental side, the degenerate Fermi gas in 
the narrow Feshbach resonant regime has been reported 
successfully in 6 Li [4l| and the Fermi-Fermi mixture of 
6 Li and 40 K j42| . Thus, it is crucially important to ex- 
plore the SO-induced superfluid physics in this regime. 
In this paper we investigate the bound state and the 
molecular signature based on the SO-driven two-channel 
model in three dimensions. Our results are given as fol- 
lows. (I) For a weak SO coupling, only one bound state 
can be found, which is identical to the result without 
SO coupling. (II) More importantly, we obtain two and 
more bound states, whose number depends strongly on 
the strength of SO coupling, in the case of the strong 
SO coupling. We also show that the bound-state num- 
ber is affected deeply by the type of SO coupling, i.e., 
the number of the ID equal Rashba-Dresselhaus SO cou- 
pling is more than that of the 2D Rashba case. (Ill) 
In addition, in the presence of SO coupling, the molec- 
ular number does not vanish (but approaches zero) even 
in the positive detuning. When varying the strength of 
SO coupling in the positive detuning, an unconventional 
molecular signature that the molecular number increases 
and then decreases is revealed. We believe that in exper- 
iments our predicted molecular behavior is a good signa- 
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ture to detect the SO-induced superfluid physics. 

Model and Hamiltonian-Yov the SO-driven two- 
channel model, the total Hamiltonian can be written for- 
mally as 

H = H F + H M + H 1 + H S . (1) 

In Hamiltonian ([I]), H-p = J2ka ^kCk^Cko- is the atom 
Hamiltonian, where C^ a is the creation operator for a 
atom with momentum k, a =t, i, £k = £k — (J-, £k = 
k 2 /(2m) is the kinetic energy, and /i is the chemical 
potential. Hm = X) q ( e q + <5o ~ 2/i)fe^6 q is the molec- 
ular Hamiltonian, where 6 q is the creation operator of 
a molecule with momentum q, e q — h 2 q 2 /(2M) with 
M = 2m is the kinetic energy of molecule, and So is 
bare detuning determined by the position S of the Fes- 
hbash resonance via a relation So = 5 + g 2 J2k l/(2 £ k). 
At lower energy, the position is given approximately by 
S ~ 2fis(B — Bo), where Bo is the magnetic field at 
which the resonance is at zero energy, and B is the tun- 
able magnetic field [43| . The atom- molecule interconver- 
sion term is governed by the following Hamiltonian H\ = 
fEqkk' ^-k+2|C k+ q t + C* 1 [, +ft C , l k , + a i 6 q , where g 
is the coupling constant that measures the amplitude of 
the decay of the molecule in the closed channel into a 
pair of the open-channel atoms. Finally, the SO coupling 
is chosen as a generalized Rashba and Dresselhaus type 
with Hamiltonian H$ = a J2 k [(k y + i\k x )C^Cui + ( k y ~ 

i\k x )C^Ck-f], where a = (an + o>d) is a generalized 
strength of SO coupling, an and a d are the strengths of 
SO coupling for Rashba and Dresselhaus types, respec- 
tively, A = (cxr — C(d)/(oir + ckd) determines the type of 
SO coupling. For example, for A = 1 (od = 0), the 2D 
Rashba SO coupling can be found. Whereas, for A = 
(ac = or), the ID equal Rashba and Dresselhaus SO 
coupling can be generated. Fortunately, this equal SO 
coupling has been realized experimentally in the ultra- 
cold neutral atoms 

In the absence of SO coupling (a = 0), Hamiltonian 
((T|) reduces to the standard two-channel model f33T - l36j ] . 
in which only the singlet Cooper paring can be formed. 
However, in the presence of SO coupling, the result is 
quite different. A crucial difference is that the singlet 
and triplet Cooper parings can coexist. As a result, the 
system exhibits rich bound states and unconventional 
molecular signature, as will be shown. For simplicity, 
we mainly focus on the case of q = in the following 
discussions. 

Two-body bound state— We begin to discuss the two- 
body bound state for the two-channel model driven by 
the generalized SO coupling from an ansatz wavefunc- 
tion of the total Hamiltonian ([1]). Due to existence of 
the singlet and triplet Cooper parings, the ansatz wave- 
function should be written formally as 

!*> = (£ /WCL^+7^)|0,0,0,0)®|0), (2) 

k aa' 



where p k > n and fa ^ (fa n and fa u ) represent the 
amplitude of probability for the singlet (triplet) Cooper 
paring, 7 represents the amplitude of probability of the 
molecule, and |0, 0,0,0) <E) |0) is the direct multiple of 
the fermion vacuum with spin flipping and molecule vac- 
uum. Substituting the wavefunction in Eq. ([3]) into the 
stationary Schrodinger equation (H —E) = 0, we find 
that the six coefficients including the energy E satisfy the 
following equations: 

37 + P k > n ak- = E k /3 Ht 

Wk'i-t + Pk't0 ak + = ^A'n , (3) 
(Pk'it + Pk'n^ ak - = E kPk'u 

{ 2{8 a -E) 1 = -g Efc'OVn + A^f) 

where = E — 2e k and k± — k y ± i\k x . Eq. ([3]) 
can not be solved directly because of lack of a coef- 
ficient equation. If we define the spin symmetry and 
anti-symmetry vectors as ip s (k) — {Pk\.-\ — Pkf-O/V^ and 
i> a {k) = (/3/4t + Pkti)/V%, the stationary Schrodinger 
equation is rewritten in the representation of tpk = 

\Pk\u Pm> PkUriPktfi 7] T as (H - E)ip k = 0, which 
leads to another equations for the coefficients f3 k > a a > 

and 7, that is, 7 = -gV2J2k' ^a(k')/(5 - E), f3 k > n = 
V2ip a (k)ak + /E k , and f3 k > u = \/2ip a (k)ak_ /Z k . Substi- 
tuting these obtained equations into Eq. ([3]) yields 

V( ^ + — ) = !tii ( 4 ) 

^ { Z 2 -4a 2 k + k_ + 2eJ g 2 ' {V 

Equation ((3J) is responsible for the bound state for 
our considered two-channel model with the generalized 
SO coupling. From Eq. ([J} we can derive the energy 
E. However, in order to prove the real bound state, 
we must cooperate with the condition E < Et, where 

E T = ^{^r- - C 1 + A )l is called the threshold en- 
ergy for the total Hamiltonian ([1} [5]. In the absence 
of SO coupling (a — 0), this threshold energy becomes 
Et = 0, as expected. In the following discussions, we 
mainly consider two interesting cases, including the 2D 
Rashba (A = 1) and ID equal Rashba and Dresselhaus 
(A = 0) SO couplings realized by the recent experiments 
[2H4| , to reveal the fundamental properties for the bound 
state. 

We first address the case of the 2D Rashba SO cou- 
pling. In the absence of SO coupling, i.e., a = 0, the 
analytical bound-state energy is derived from Eq. Q by 
E = (- 9 4 m 3 + 32tt 2 S - g 2 y/g 2 m e - 64m 3 7r 2 J)/(327r 2 ). 
It implies that in such a case only one bound state can 
be found [32j]. In the presence of SO coupling, the ex- 
plicit expression for the bound-state energy can not be 
obtained. However, for a weak SO coupling, Eq. @ 
is simplified as 2mi{E + ma 2 ) / (8ir^/^E) = (E-5)/g 2 
with the help of a Taylor expansion with respect to the 
strength of SO coupling. In this case, Hamiltonian ([1]) 
also exhibits one bound state, as shown in Fig. 1(a). 
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FIG. 1: (Color online) The bound-state energy for the 2D 
Rashba SO coupling (A = 1) as a function of the detuning 5 
for the different strengths of SO coupling (a) otK F = l.OEp, 
(b) aK F = 3.6£f and (c) aK F = 7.0E F when T = 0.1. In Fig. 
1(a), the red open symbol corresponds to the analytical result 
(AR) and the black solid line represents the direct numerical 
simulation (NS). 



By further solving the above equation approximately, we 
find that the bound-state energy is proportional to —a 4 . 
This behavior agrees well with the numerical simulation, 
as also shown in Fig. 1(a), and implies that the bound- 
state energy can decrease by increasing the strength of 
SO coupling. For the strong SO coupling, the perturba- 
tion method is invalid. Here we numerically solve Eq. (|4]) 
to evaluate the energy E. Quite surprisingly, Hamilto- 
nian |l| has two bound states in the strong SO coupling, 
as shown in Fig. 1(b). Moreover, the bound-state num- 
ber depends strongly on the strength of SO coupling , i.e., 
when increasing the strength of SO coupling, the bound- 
state number increases, as shown in Fig. 1(c). The ex- 
istence of the multi-bound states is a new phenomenon 
induced by SO coupling and may generate exotic super- 
fluid physics. 

Another important observation of Eq. (0| is that the 
bound-state number is affected deeply by the type of SO 
coupling. In Fig. 2, we plot the bound-state energy 
with respect to the detuning 5 for the different types 
of SO coupling including (a) A = 1 (2D Rashba) and (b) 
A = (ID equal Rashba and Dresselhaus). It can be seen 
clearly that for a fixed strong SO coupling, the bound- 
state number for the ID equal Rashba and Dresselhaus 
SO coupling is more than that of the 2D Rashba SO 
coupling. In order to prove this argument again, we nu- 
merically calculate the bound-state number in the range 
of the strength aKp of SO coupling from to 6AEp 
when r = 0.1 and 6 = —3Ep. Wc find approximately 
that for the 2D Rashba SO coupling, the bound-state 
number is given by N R = 1 for aKp/Ep £ [0.0,3.0] and 
N R = 2 for aKp/Ep G [3.1,6.4]. Whereas, for the ID 
equal Rashba and Dresselhaus SO coupling, the bound- 
state number becomes Nrd = 1 for aKp/Ep € [0.0, 2.9], 



FIG. 2: (Color online) The bound-state energy as a function of 
the detuning S for the different types of SO coupling including 
(a) A = 1 (2D Rashba) and (b) A = (ID equal Rashba and 
Dresselhaus), when F = 0.1 and aK F = 4.2E F . 



ATrd = 2 for aKp/Ep G [3.0,3.5], AT RD = 3 for 
aKp/Ep £ [3.6,3.9], N RD = 4 for aKp/Ep G [4.0,4.9], 
A/rx) = 5 for aKp/Ep £ [5.0,5.4], and N RU = 6 for 
aKp/Ep £ [5.5,6.4]. 

Molecular signature-After discussing the bound state 
driven by SO coupling in the two-channel model, now we 
explore the experimentally-measurable molecular signa- 
ture. In terms of the Hellmann-Feymann theorem, the 
molecular number is obtained by 

N = (*\b%\*) = (*\ — \*) = —, (5) 

where the bound-state energy E can be derived from Eq. 
0, together with the condition E < Et- Here we are 
interested in the regime for the one bound state. 

In Fig. 3 we plot the scaled molecular number 2N$/N 
of both the 2D Rashba and the ID equal Rashba and 
Dresselhaus SO couplings with respect to the detuning 5 
for the different strengths of SO coupling. In the absence 
of SO coupling (aKp = 0), the striking molecular signa- 
ture can be found in the case of the negative detuning 
(S < 0). However, for the positive detuning (5 > 0), the 
physical bound state disappears 44], and thus there is 
no real molecular number, as shown in black solid line 
of Fig. 3. With the increasing of the strength of SO 
coupling, the molecular number in the negative detun- 
ing decreases, as shown in red dash line of Fig. 3. The 
physics can be understood as follows. In the two chan- 
nel model, the molecules play two roles. One is that 
they interact directly with the atoms via Hamiltonian 
Hj. The other (the most important) is that they induce 
the indirect atom-atom interactions, which generate the 
Cooper pairing. When the SO coupling enhances the 
Cooper pairing the molecules are thus suppressed 

because the system need guarantee a conserved number 

N = 2blb + J^ka c ia C *°- 0n the other hand > in the 
presence of SO coupling, there still exists a large molec- 
ular signature in the positive detuning, which is contrast 
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FIG. 3: (Color online) The molecular number of both (a) 
the 2D Rashba SO coupling with A = 1 and (b) the ID equal 
Rashba and Dresselhaus SO coupling with A = as a function 
of the detuning 5 for the different strengths of SO coupling 
cxKf = 0.0i?F (Black solid line), cxKf — 1.5E? (Red dashed 
line) and aKp — 3.0-Ef (Cyan dash-dotted line) when V — 
0.1. 
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FIG. 4: (Color online) The molecular number of both (a) 
the 2D Rashba SO coupling with A = 1 and (b) the ID equal 
Rashba and Dresselhaus SO coupling with A = as a function 
of the strength aKp of SO coupling for the different detung- 
ings S — S.OEp (Black solid line), 5 — 1.5-Ep (Red dashed 
line) and 8 — 3.0-Ep (Cyan dash-dotted line) when T = 0.1. 



to the result without SO coupling. Finally, by means of 
the perturbation method, the molecular number is ob- 
tained as 2N /N = C\+ C*2a 2 , where G\ and C2 are 
the complicate coefficients independent of the strength 
of SO coupling. In the absence of SO coupling, C\ = 



for 5 > and C\ ~ 1 for 8 < 0. This analytical expression 
shows that the molecular number does not vanish (but 
approaches zero) for a positive detuning in the presence 
of SO coupling, as shown in cyan dash-dotted line of Fig. 
3. 

In Fig. 4, we plot the molecular number of both the 
2D Rashba and the ID equal Rashba and Dresselhaus 
SO couplings with respect to the strength of SO cou- 
pling for the different detunings. In the negative detun- 
ing (8 < 0) we shows again that the SO coupling sup- 
presses the molecular number (black solid line of Fig. 4). 
However, for the zero and positive detunings (8 > 0), we 
find that the molecular number increases and then de- 
creases, and eventually, approaches zero, when increas- 
ing the strength of SO coupling. Moreover, with the 
increasing of the negative detuning, the corresponding 
peak of the molecular number decreases. In contrast to 
the pairing gap that varies smoothly as a function of 
the strength of SO coupling [6j-|8[ , we believe that in ex- 
periments our predicted unconventional behavior for the 
molecular number is a good signature to detect the SO- 
induced superfiuid physics. 

Conclusions-hi summary, motivated by the recent ex- 
perimental developments about the degenerate Fermi 
gas, we have investigated the SO-driven BCS-BEC 
crossover in the narrow Feshbash resonance based on the 
generalized two-channel model. We have found that the 
SO coupling can generate rich bound states and uncon- 
ventional molecular signature. In experiments, we believe 
that our predicted molecular behavior is a good signature 
to detect the SO-induced superfiuid physics. In addition, 
since the multi-bound states can promote fluctuations 
and possibly inhibit condensation, some SO-induced ex- 
otic phenomena may occur in the narrow Feshbash reso- 
nance. 
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Note added-During preparing this paper, we noticed 
that two bound states for the SO-driven two-channel 
model was also predicted by V. B. Shenoy in terms of 
a renormalizable quantum field theory j45| . Our results 
can agree with each other where they overlap. 
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